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IMPAIRMENT CORRELATION ESTIMATION 
IN A SPREAD SPECTRUM SYSTEM 



BACKGROUND OF THE INVENTION 
5 The present invention relates generally to signal processing parameters for 

processing received signals in a wireless communication system, and more 
particularly to estimating signal impairment correlations of multiple images of the 
received signals. 

Wireless receivers in spread spectrum systems, such as code division 

10 multiple access (CDMA) systems, are well known in the communication arts. 

Typically, wireless receivers operating in a spread spectrum system include a RAKE 
receiver, which is responsible for demodulating multiple signal images of a received 
multi-path signal. In general, a RAKE receiver in a base station or a mobile station 
receives a desired signal in the presence of other interfering signals intended for or 

1 5 sent from other base stations and/or mobile stations. By using calculated weighting 
factors to coherently combine correlated symbols corresponding to the desired 
received signal, the RAKE receiver uses multi-path reception to improve the signal- 
to-noise ratio (SNR) of the received multi-path signal. 

While conventional RAKE receivers are effective at reducing channel path 

20 noise, receiver noise, etc., and therefore, are effective at improving the SNR of 
received multi-path signals, conventional RAKE receivers are typically ineffective 
against interference noise caused by self-interference, multi-user access 
interference, etc. To address noise and interference, a generalized RAKE (G-RAKE) 
receiver may be used instead of the conventional RAKE receiver, as discussed in co- 

25 assigned U.S. Patent No. 6,363,104 to Bottomley entitled "Method and Apparatus for 
Interference Cancellation in a RAKE Receiver," issued 26 March 2002. This patent is 
incorporated herein by reference. Like a RAKE receiver, a G-RAKE receiver 
correlates received signal images, weights the correlated signal images using 
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weighting factors, and coherently combines the weighted signal images. However, 
unlike the RAKE receiver, the weighting factors of the G-RAKE are based on 
estimated impairment correlations. To suppress noise and interference, the G-RAKE 
combiner uses impairment correlation based weighting factors to weight the 
5 correlated signal images and to coherently combine the weighted images. As 

discussed further below, the success of the interference suppression in a G-RAKE 
receiver is dependent on the accurate estimation of these impairment correlations. 

SUMMARY OF THE INVENTION 

10 The present invention comprises a method and apparatus for deriving an 

impairment correlation matrix to process signals received at a wireless receiver over 
multiple paths of a multi-path channel. The receiver despreads symbols of the 
received signals and determines channel estimates from the despread symbols. 
First and second impairment correlation estimators estimate first and second 

15 impairment correlation matrices based on the despread symbols. Based on the first 
and second impairment correlation matrices, a correlation processor derives the 
impairment correlation matrix. 

The receiver may use the derived impairment correlation matrix to generate 
one or more processing parameters for processing the received signals. In one 

20 embodiment, a weight calculator calculates weighting factors based on the derived 
impairment correlation matrix. The calculated weighting factors are then used to 
coherently combine the despread symbols to suppress noise and interference in the 
received signal. In another embodiment, a signal-to-interference ratio (SIR) 
calculator calculates an SIR estimate based on the derived impairment correlation 

25 matrix. The estimated SIR may be used for rate adaptation or as part of the transmit 
power control for the base station or mobile station. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 illustrates an exemplary multi-path wireless communication system. 
FIG. 2 illustrates an exemplary wireless communication system receiver. 
FIG. 3 illustrates an exemplary baseband processor of the present invention. 
5 FIG. 4 illustrates an exemplary parameter processor of the present invention. 

FIG. 5A illustrates an exemplary impairment correlation estimator according 
to the present invention. 

FIG. 5B illustrates another exemplary impairment correlation estimator 
according to the present invention. 
10 FIG. 6 illustrates an exemplary parametric estimator. 

FIG. 7 illustrates an exemplary non-parametric estimator. 
FIGs. 8A - 8D illustrate exemplary correlation processors for the impairment 
correlation estimators of Figures 5A and 5B. 

FIG. 9 illustrates an exemplary method for deriving an impairment correlation 
1 5 matrix according to the present invention. 

FIG. 10 illustrates another exemplary impairment correlation estimator 
according to the present invention. 

FIG. 11 illustrates an exemplary correlation processor for the impairment 
correlation estimator of Figure 10. 
20 FIG. 12 illustrates another exemplary impairment correlation estimator 

according to the present invention. 

FIG. 13 illustrates an exemplary correlation processor for the impairment 
correlation estimator of Figure 12. 

FIG. 14 illustrates another exemplary method for deriving an impairment 
25 correlation matrix according to the present invention. 
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DETAILED DESCRIPTION OF THE INVENTION 
Figure 1 illustrates an exemplary spread spectrum wireless communication 
system 10. The wireless communication system 10 includes at least one base 
station 20, at least one mobile station 30, and possibly one or more interfering 
5 objects 40. As used herein, the term "mobile station" may include a cellular 
radiotelephone with or without a multi-line display; a Personal Communications 
System (PCS) terminal that may combine a cellular radiotelephone with data 
processing, facsimile and data communications capabilities; a personal data 
assistant (PDA) that can include a radiotelephone, pager, Internet/intranet access, 

10 Web browser, organizer, calendar, and/or a global positioning system (GPS) 

receiver; and a conventional laptop and/or palmtop receiver or other appliance that 
includes a radiotelephone transceiver. Mobile stations may also be referred to as 
"pervasive computing" devices. 

Base station 20 includes one or more antennas 22 for transmitting/receiving 

15 spread spectrum signals with one or more symbols to/from mobile station 30. The 
transmitted signals typically include traffic and pilot signals. Objects, such as 
interfering object 40, cause multiple echoes of the transmitted symbols to arrive at 
mobile station 30 at different times. Receiver 34 processes the multiple signal 
images at mobile station 30. Similarly, mobile station 30 may transmit symbols via 

20 one or more antennas 32 along multiple paths to base station 20, where receiver 24 
processes the multiple received symbol images. 

Figure 2 illustrates an exemplary receiver 24, 34 for base station 20 and/or 
mobile station 30. The receiver 24, 34 receives and processes the symbols of the 
received signals to generate received symbol estimates. An exemplary receiver 24, 

25 34 includes a front end 50, a baseband processor 100, and an additional processor 
52. Front end 50 typically includes filters and/or mixers to convert the received 
signals to baseband signals, and conversion circuits, such as analog-to-digital 
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converters, to produce a series of digitized baseband signal samples r(t) 
corresponding to the baseband signal. Baseband processor 100 demodulates the 
baseband signal samples r(t) to produce symbol estimates corresponding to the 
received signal. These symbol estimates are then processed further, as necessary, 
5 in additional processor 52. For example, additional processor 52 may include a turbo 
decoder (not shown) that determines information bit values based on the symbol 
estimates provided by the baseband processor 100. These information bit values 
may then be converted to speech, images, etc. 

Baseband processor 100 may include a G-RAKE receiver, as shown in Figure 

10 3, for determining the desired symbol estimates while suppressing noise and 

interference. An exemplary G-RAKE receiver includes one or more correlation units 
102, a combining circuit 104, a finger placement processor 106, and a parameter 
processor 110. Based on the baseband signal r(t), finger placement processor 106 
identifies and selects one or more delays 6 k corresponding to signal paths of the 

15 multi-path channel as well as extra interference suppression delays. Each selected 
delay d* is provided to correlation unit 102, which despreads one or more traffic 
channel signals present in the baseband signal r(t) based on the provided delays 6 k . 
As used herein, traffic channel signals include broadcast, control, and/or individual 
data signals. The selected delays d k are also provided to parameter processor 110 

20 for computing a weighting factor w* for each delay d* and, optionally, for computing 
an SIR estimate. Combining circuit 104 receives each despread traffic channel 
signal y k and the corresponding weighting factor w*, and generates symbol estimates 
s(m) by weighting and coherently combining the despread traffic channel signals y k 
according to methods known in the art. 

25 As shown in Figure 3, parameter processor 110 calculates weighting factors 

Wfr based on the baseband signal r(t) and on the selected delays d*. Figure 4 
illustrates an exemplary parameter processor 1 10 according to the present invention. 
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Parameter processor 110 includes a pilot correlation unit 1 12, a modulation remover 
1 14, a channel tracker 1 16, a weight calculator 1 18, an impairment correlation 
estimator 120, and optionally, an SIR calculator 119. The pilot correlation unit 112 
uses the delays d* provided by the finger placement processor 106 to despread one 
5 or more pilot channel signals present in the baseband signal r(t). Modulation 

remover 114 removes the symbol modulation from the despread pilot symbols p* by, 
e.g., multiplying each despread pilot symbol p k by the complex conjugate of the 
known value to provide demodulated and despread pilot symbols. In systems with 
higher order modulation, modulation remover 114 may also divide the despread pilot 

10 symbol p k by the pilot symbol value, as is known in the art. Data symbols and 
decision feedback may also be used. 

Channel tracker 116 uses the demodulated and despread pilot symbols to 
estimate the channel response coefficients, and therefore, to produce channel 
estimates. Based on the demodulated and despread pilot symbols from modulation 

15 remover 114 and the channel estimates from channel tracker 116, impairment 

correlation estimator 120 estimates an impairment correlation matrix R having one or 
more impairment correlation values corresponding to the set of delay paths being 
processed by receiver 24, 34. Note that the impairment correlation values in 
impairment correlation matrix R are based on both noise correlation terms and 

20 interference correlation terms. Because those skilled in the art will appreciate that 
"covariance" is a special case of "cross-correlation" with zero mean, the terms 
"correlation" and "covariance," as used herein, should be understood as 
interchangeable unless the context of a particular passage makes an explicit 
distinction between the two terms. 

25 Weight calculator 118 uses the channel estimates and the impairment 

correlation values in the impairment correlation matrix R to calculate weighting 
factors corresponding to the delay paths of receiver 24, 34. The weighting factors 
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are used to combine traffic despread values to form symbol estimates. As discussed 
above, accurate weighting factors w* are important for accurate interference 
suppression in G-RAKE receivers. Because the weighting factors w* directly relate to 
the estimated impairment correlation matrix R, interference suppression in G-RAKE 
5 receivers depends on the accurate estimation of the impairment correlation matrix R. 

Figure 5A illustrates an exemplary impairment correlation estimator 120 for 
estimating impairment correlation matrices according to the present invention. 
Impairment correlation estimator 120 comprises a first impairment correlation 
estimator 130, a second impairment correlation estimator 140, and a correlation 

10 processor 160. First impairment correlation estimator 130 generates a first 

impairment correlation matrix R A based on the channel estimates and despread 
values. Similarly, the second impairment correlation estimator 140 generates a 
second impairment correlation matrix R B based on the channel estimates and 
despread values. Correlation processor 160 then derives the final impairment 

15 correlation matrix R based on the first and second impairment correlation matrices, 
R A and R B . 

The first and second impairment correlation estimators 130, 140 may be any 
impairment estimator that estimates an impairment correlation matrix. Exemplary 
first and second impairment correlation estimators 130, 140 may include those 

20 described in U.S. Patent No. 6,363,104 ("Method and Apparatus for Interference 
Cancellation in a Rake Receiver'' issued 26 March 2002), U.S. Patent Application 
Serial No. 09/344,899 ("RAKE Combining Methods and Apparatus Using Weighting 
Factors Derived from Knowledge of Spread Spectrum Signal Characteristics" filed 25 
June 1999), and U.S. Patent Application Serial No. **/***,*** ("Method and Apparatus 

25 for Parameter Estimation in a Generalized RAKE Receiver", filed 12 March 2004, 
attorney docket number 4015-5161), all of which are herein incorporated by 
reference. Additional exemplary impairment correlation estimators are also 
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described in "Approaches for Fast Adaptive Generalized RAKE Reception" to 
Bottomley et al., in Research Disclosure, November 2003; "Performance of CDMA 
Mobile Communication Systems Using Antenna Arrays" to Suard et al., in Proc. IEEE 
ICASSP, 27-30 April 1993, pp. IV-153 through IV-156; and "Pilot Channel-aided 
5 Techniques to Compute the Beamforming Vector for CDMA Systems with Antenna 
Array" to Choi in IEEE Trans. Veh. Techno!., vol. 49, pp. 1760-1775, September 
2000. 

In one exemplary embodiment of the present invention, the first and second 
impairment correlation estimators 130, 140 may be equivalent impairment correlation 

10 estimators programmed with different settings. Alternatively, in another exemplary 
embodiment of the present invention, the first and second impairment correlation 
estimators 130, 140 may be different types of impairment correlation estimators. For 
example, as shown in Figure 5B, the first impairment correlation estimator 130 may 
be a parametric estimator 130 and the second impairment correlation estimator 140 

15 may be a non-parametric estimator 140. These estimators share an initial 

processing step in which channel estimates and modulation-removed despread 
values are provided to signal remover 122. Signal remover 122 subtracts the 
channel estimates from the modulation-removed despread values to form an error 
vector e. In this embodiment, the correlation processor 160 derives the final 

20 impairment correlation matrix R from a parametric impairment correlation matrix R A 
and a non-parametric impairment correlation matrix R B . While the following 
discussions describe the impairment correlation estimator 120 of the present 
invention in terms of a parametric estimator 130 and a non-parametric estimator 140, 
it will be appreciated by those skilled in the art that the first and second impairment 

25 correlation estimators 130, 140 may comprise any impairment correlation estimators 
known in the art. 
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To facilitate further discussion of the present invention, Figures 6 and 7 
illustrate exemplary parametric and non-parametric estimators 130, 140, respectively. 
A parametric estimator 130 models one or more sources of interference based on the 
received signal and generates the parametric impairment correlation matrix R A based 
5 on these interference models. As described in U.S. Patent Application Serial No. 
**/***,*** ("Method and Apparatus for Parameter Estimation in a Generalized RAKE 
Receiver", filed 12 March 2004, attorney docket number 4015-5161) and shown in 
Figure 6, an exemplary parametric estimator 130 comprises a correlation computer 
132, a structure element computer 134, a scaling parameter estimator 136, and an 

10 impairment correlation calculator 138. Correlation computer 132 receives the error 
vector e from the signal remover 122 and generates correlation measurements 
based on the values in the error vector e. In general, correlation computer 1 32 
multiplies each error value in error vector e by the conjugates of another error value 
to produce the correlation measurements. These correlation measurements may, for 

1 5 example, be the result of several products averaged together over a CDMA time slot. 

Structure element computer 134 constructs "structure" elements that define 
the structure of the parametric impairment correlation model(s). A typical structure 
element computer 134 constructs an interference structure element R, and a noise 
structure element R 2 . The structure elements R^ and R 2 are calculated according to 

20 the parametric impairment model(s) by structure element computer 134 based on the 
channel estimates provided by the channel tracker 116. More specifically, the values 
of the interference structure element R^ are based on the channel estimates and 
knowledge of the receiver pulse shape, while the noise structure element R 2 is made 
of values that depend on the known receiver pulse shape. 

25 Based on the correlation measurements from correlation computer 132 and 

the structure elements R 1a R 2 from the structure element computer 134, scaling 
parameter estimator 136 derives scaling parameters, a and /?, for each of the 
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structure elements. Using the scaling parameters a and fi, and structure elements 
R t and R 2 , the impairment correlation calculator 138 combines the structure elements 
to estimate the parametric impairment correlation matrix R A . According to an 
exemplary parametric estimator 130, impairment correlation calculator 138 combines 
5 the structure elements according to: 

R A =aR 1 + pR 2 . (Eq. 1) 

Because the parametric impairment correlation matrix R A depends on the 
channel estimates, parametric estimator 130 is able to track rapid variations in the 
modeled interference, particularly in own-cell interference noise correlation. 

10 However, because the parametric estimator 130 does not model all forms of 

interference, some impairment correlations may not be properly accounted for in the 
impairment correlation matrix R A . 

To address this issue, the second impairment correlation estimator 140 of 
Figure 5B may comprise a non-parametric estimator 140. Non-parametric estimator 

15 140, as shown in Figure 7, includes a correlation computer 142 and a smoothing filter 
144. As with the parametric estimator 130, the correlation computer 142 in the non- 
parametric estimator 140 generates correlation measurements based on the error 
values in error vector e. The correlation measurements are then filtered over time, 
e.g., over multiple CDMA time slots of the received CDMA signal, in smoothing filter 

20 144 to produce the non-parametric impairment correlation matrix R B . 

Because the accuracy of the non-parametric impairment correlation matrix R B 
is directly related to the length of time used to smooth the correlation measurements, 
the non-parametric estimator 140 cannot accurately track rapid variations in the 
interference. However, because the non-parametric impairment correlation matrix R B 

25 is not dependent on particular interference models, and instead is directly derived 

from the received signal, it necessarily includes all forms of interference. As such, by 
deriving an impairment correlation matrix R based on the parametric and the non- 



10 of 37 



Ericsson Ref. No.: P18409-US1 
C&B Ref. No.: 4015-5177 



parametric impairment correlation matrices, R A and R B , the impairment correlation 
estimator 120 of the present invention produces an impairment correlation matrix R 
that can track rapid variations in the interference while simultaneously including all 
forms of interference present in the received signal. 
5 Correlation processor 160 may derive the final impairment correlation matrix 

based on the parametric and non-parametric impairment correlation matrices 
according to any preferred method. In a first embodiment, correlation processor 160 
may comprise a selector 160, as shown in Figure 8A, that selects either the 
parametric impairment correlation matrix R A or the non-parametric impairment 
10 correlation matrix R B for each CDMA time slot dependent upon which impairment 
correlation matrix best satisfies a predetermined criteria, such as interference 
suppression. For example, an SIR metric, M(R), may be calculated for each 
impairment correlation matrix, R A and R B , according to: 

MfRJ^R-'c, (Eq. 2) 

15 where c is a vector of the channel estimates and where c H represents the Hermitian 
transpose of the vector of channel estimates. Selector 160 may select the 
impairment correlation matrix that maximizes the SIR metric M(R). 

In another embodiment, shown in Figure 8B, the correlation processor 160 
may comprise a smoothing filter 162, a subtractor 164, and a summer 166. In this 

20 embodiment, smoothing filter 162 filters the values in the parametric impairment 
correlation matrix R A to produce a filtered parametric impairment correlation matrix 
R A that has a filtering similar to that used by the non-parametric impairment 
correlation estimator 140. Subtractor 164 subtracts the filtered parametric 
impairment correlation matrix R A from the non-parametric impairment correlation 

25 matrix Rb to generate a residual correlation matrix Rc. The residual correlation 
matrix Rc comprises impairment correlation values that correspond to impairment 
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correlations detected by the non-parametric estimator 140, but that are not modeled 
by the parametric estimator 130. By adding the residual correlation matrix R c to the 
parametric impairment correlation matrix R A in summer 166, a final impairment 
correlation matrix R that includes both modeled and non-modeled interference is 
5 derived. 

If the parametric estimator 130 does not model a significant interference 
source, the un-modeled interference source may bias the estimation of the scaling 
parameters, a and /?. Correlation processor 160 may mitigate this bias effect by 
feeding back a term for the un-modeled interference source to the parametric 
10 estimator 130 and improving the interference model based on this term. To 

implement this mitigation effect, correlation processors 160 of Figures 8C and 8D 

include an additional smoothing filter 168 for generating a filtered or final version R c 
of the residual correlation matrix R c . By feeding the filtered residual correlation 
matrix R c back into the parametric estimator 130, the filtered residual correlation 

15 matrix R c augments the parametric impairment correlation matrix R A . For example, 
the filtered residual correlation matrix R c may be provided to both the impairment 
correlation calculator 138 and the scaling parameter estimator 136 of the parametric 
estimator 130. In response, the scaling parameter estimator 136 may estimate a and 
fi t in addition to an additional scaling factor y The impairment correlation calculator 

20 138 then augments the parametric impairment correlation matrix R A according to: 

R' A = aR, + pR 2 + yR c , (Eq. 3) 

where R' A represents the augmented parametric impairment correlation matrix. 
Note that parametric estimator 130 uses the augmented matrix R' A of Equation 3 to 
determine a, /?, and y These fitting parameters are then used in Equations 1 and 3 

25 to determine R A and R' A . For the correlation processor 160 of Figure 8C, R A , 
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calculated in Equation 1 and based on the new a and J3 provided by scaling 
parameter estimator 136, is summed with residual correlation matrix R c in summer 
166 to generate the final impairment correlation matrix R. For the correlation 
processor 160 of Figure 8D, the augmented matrix R' A may be used as the final R 
5 matrix. 

Figure 9 illustrates an exemplary method for implementing the impairment 
correlation estimator 120 of Figure 5B using the correlation processor 160 of Figures 
8C or 8D. The impairment correlation estimation process begins (block 205) by 

initializing the filtered residual correlation matrix R c to zero (block 210). The 
10 parametric estimator 130 estimates the parametric impairment correlation matrix R A 
and the "augmented" parametric impairment correlation matrix R' A as described 
above using Equation 1 and Equation 3 (block 215), and the non-parametric 
estimator 140 estimates the non-parametric impairment correlation matrix R B as 
described above with reference to Figure 7 (block 220). Because the filtered residual 
15 correlation matrix R c has been initialized to zero, the initial estimation of the 
augmented impairment parametric correlation matrix R' A is the same as the 
parametric impairment correlation matrix R A of Equation 1. 

After filtering the parametric impairment correlation matrix R A (block 225), 
the residual correlation matrix R c is generated by subtracting the filtered parametric 
20 impairment correlation matrix R A from the non-parametric impairment correlation 

matrix R B (block 230). Smoothing filter 168 then filters the residual correlation matrix 
Rc (block 235) to generate the filtered residual correlation matrix R c that is fed back 
to the parametric estimator 130 and that replaces the initialized filtered residual 
correlation matrix R c . 
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The steps of generating the augmented and non-augmented parametric 
impairment correlation matrices R' A and R A (block 215), the non-parametric 
impairment correlation matrix R B (block 220), the residual correlation matrix Rc 
(blocks 225 and 230), and the filtered residual correlation matrix R c (block 235) are 
5 repeated until a stopping criteria has been met (block 240). For example, the 
stopping criteria may simply be defined by a desired number of iterations, such as 
two iterations or three iterations. Alternatively, the stopping criteria may be based on 

an evaluation of the filtered residual correlation matrix R c to see if R c is converging 
and/or on an evaluation of the scaling factors or, /?, and ^to see if they are 

10 converging. In any event, once a sufficient number of iterations have been 

performed (block 240), the final impairment correlation matrix R is calculated as the 
augmented parametric impairment correlation matrix R' A (block 245), as shown in 
Figure 8D, and the process ends (block 250). This process is repeated for one or 
more time slots of the received signal. Note that alternatively, the final impairment 

15 correlation matrix R may be calculated by adding the parametric impairment 

correlation matrix R A to the residual correlation matrix R c , as is done in Figure 8C. 

In another exemplary embodiment of the present invention, the impairment 
correlation estimator 120 may be implemented as shown in Figure 10. The 
impairment correlation estimator 120 of Figure 10 includes the signal remover 122, 

20 the parametric estimator 130, a non-parametric estimator system 150, correlation 
processor 160, and optionally, a control unit 124. As described above, parametric 
estimator 130 generates a parametric impairment correlation matrix R A from the error 
vector e provided by the signal remover 122. However, unlike the non-parametric 
estimator 140 described above, the non-parametric estimator system 150 of Figure 

25 10 generates a non-parametric impairment correlation matrix Ro based on the 

parametric impairment correlation matrix R A . Correlation processor 160 then derives 
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the final impairment correlation matrix from the parametric impairment correlation 
matrix R A and the non-parametric correlation matrix R D generated by the parametric 
estimator 130 and the non-parametric estimator system 150, respectively. 

As shown in Figure 10, non-parametric estimator system 150 includes 
5 whitening matrix unit 152, a matrix multiplier 154, and a non-parametric estimator 
140. Whitening matrix unit 152 computes a whitening matrix F" 1 from the parametric 
impairment correlation matrix R A . In exemplary embodiments, the whitening matrix 
F~ 1 is derived from the inverse of the square root of the parametric impairment 
correlation matrix R A , where the square root is obtained through Cholesky 
10 factorization as well understood in the art. Specifically, if the parametric impairment 
correlation matrix R A is factored using the Cholesky decomposition, R A may be 
represented as: 

Ra=FF h , (Eq.4) 
where the "H" represents the Hermitian transpose operation. The whitening matrix is 

15 the inverse of F, which is denoted as F 1 . Techniques for obtaining a whitening 

matrix are known in the art, and are, for example, described further in section 6.6 of 
"Detection of Signals in Noise, 2 nd Edition" by R.N. McDonough and A.D. Whalen, 
published by Academic Press in 1995. It will be appreciated that when R A is 
diagonally dominant, a simple form of Cholesky factorization may be used, which 

20 takes the square root of each diagonal element to generate the matrix F. The 
whitening matrix F" 1 is then obtained by taking the reciprocal of each diagonal 
element. Whitening matrix unit 152 may determine if the parametric impairment 
correlation matrix R A is diagonally dominant according to any means known in the 
art. For example, whitening matrix unit 152 may determine that parametric 

25 impairment correlation matrix R A is diagonally dominant by comparing the sum of the 
magnitude square of the off-diagonal elements to a threshold. If the sum exceeds 
the threshold, R A is not diagonally dominant. Alternatively, the whitening matrix unit 
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152 may compute a first SIR using R A with the diagonal elements and a second SIR 
using the R A without the diagonal elements. Comparing the ratio of the first and 
second SIRs to a threshold determines if the parametric impairment correlation 
matrix R A is diagonal dominant. 
5 In any event, matrix multiplier 154 multiplies the error vector e by the 

whitening matrix F' 1 according to 

x = F 1 e (Eq. 5) 

to produce a whitened error vector x. Non-parametric estimator 140, which 
corresponds to the non-parametric estimator 140 of Figure 7, then estimates the non- 
10 parametric impairment correlation matrix R D based on the whitened error vector x. 

Correlation processor 160 derives the final impairment correlation matrix R 
based on the parametric impairment correlation matrix R A , the whitening inverse 
matrix F, and the non-parametric impairment correlation matrix R D . To achieve this 
derivation, an exemplary correlation processor 160 may comprise a converter 170 
15 and a combiner 172, as shown in Figure 11. Converter 170 applies the whitening 
matrix inverse F to the non-parametric impairment correlation matrix R D to produce 
the final residual correlation matrix R c . In exemplary embodiments, the converter 
170 generates the final residual correlation matrix R c according to: 

Rc =F(R D -l)F H . (Eq.6) 
20 Combiner 172 combines the parametric impairment correlation matrix R A with the 
final residual correlation matrix R c by summing the two matrices to derive the final 
impairment correlation matrix R. 

Alternatively, as shown in Figure 12, correlation processor 160 may derive the 
final impairment correlation matrix R based on an augmented parametric impairment 
25 correlation matrix R' A , the whitening matrix inverse F, and the non-parametric 

impairment correlation matrix Ro- In this embodiment, the correlation processor 160 
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of Figure 12 comprises a converter 170, as shown in Figure 13. As with the 
converter 170 of Figure 1 1 , converter 170 of Figure 13 applies the whitening matrix 
inverse F to the non-parametric impairment correlation matrix R D to generate a final 

residual correlation matrix R c . Final residual correlation matrix R c is then fed back 

5 into the parametric estimator 1 30. Parametric estimator 1 30 uses the final residual 

correlation matrix R c to generate an augmented parametric impairment correlation 

matrix R' A , which may be derived from the final residual correlation matrix R c 
according to Equation 3, as described above. Correlation processor 160 then 
derives the final impairment correlation matrix R directly from the augmented 
10 parametric impairment correlation matrix R' A . 

Figure 14 illustrates an exemplary method for implementing the impairment 
correlation estimator 120 of Figure 12 using the correlation processor 160 of Figure 
13. The impairment correlation estimation process begins (block 305) by initializing 

the final residual correlation matrix R c to zero (block 310). Parametric estimator 130 
15 forms the parametric impairment correlation matrix R A and the "augmented" 

parametric impairment correlation matrix R' A by learning fitting parameters and 
applying these fitting parameters as described above using Equations 1 and 3, 
respectively, (block 315). Whitening matrix unit 152 generates the whitening matrix 
F" 1 and whitening matrix inverse F from the parametric impairment correlation matrix 
20 R A as described above (block 320). Matrix multiplier 154 and the non-parametric 
estimator 140 generate the non-parametric impairment correlation matrix Rq as 
described above (block 325). Converter 170 applies the whitening matrix inverse F 
to the non-parametric impairment correlation matrix R D to generate the final residual 

correlation matrix R c (block 330), replaces the initialized final residual correlation 
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matrix R c with the newly calculated final residual correlation matrix R c (block 335), 
and provides R c to the parametric estimator 130. 

The steps of generating the parametric impairment correlation matrices R A 
and R' A (block 315), the non-parametric impairment correlation matrix R D (blocks 

5 320 and 325), and the final residual correlation matrix R c (blocks 330 and 335) are 
repeated until a stopping criteria has been met (block 340). For example, the 
stopping criteria may simply be defined by a desired number of iterations, such as 
two iterations or three iterations. Alternatively, the stopping criteria may be based on 
an evaluation of the final residual correlation matrix R c to see if R c is converging 
10 and/or on an evaluation of the scaling factors a, /?, and /to see if they are 
converging. In any event, once a sufficient number of iterations have been 
performed (block 340), the final impairment correlation matrix R is defined as the 
augmented parametric impairment correlation matrix R' A (block 345), and the 

process ends (block 350). This process is repeated for one or more time slots of the 

15 received signal. 

As shown in Figures 10 and 12, some embodiments of the impairment 
correlation estimator 120 may include a control unit 124. Control unit 124 may 
simplify the whitening matrix calculations, as discussed above. In addition, control 
unit 124 may conserve power associated with the base station 20 or mobile station 

20 30 by intelligently applying the impairment correlation matrix estimation process of 
the present invention only to the time slots where the parametric impairment 
correlation matrix estimate does not satisfy predetermined color criteria. For 
example, for each time slot of a received spread spectrum signal, control unit 124 
evaluates the current parametric impairment correlation matrix R A (non-augmented) 

25 to determine how much color is present in the current impairment correlation matrix 
R A of the current time slot. When control unit 124 determines that the current 
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impairment correlation matrix R A of a particular time slot contains sufficient color, as 
compared to a predetermined color criteria, control unit 124 temporarily deactivates 
non-parametric estimator system 150 and the converter 170 in correlation processor 
160 to conserve power, and identifies the parametric impairment correlation matrix 
5 R A as the final impairment correlation matrix R for that time slot. One exemplary 
color criteria may be a fitness measure between R A and the measured correlations 
that is used to decide if further sources of correlation may be learned by system 150. 
For example, if the fitness measure indicates that R A sufficiently models the 
measured correlations, then the non-parametric estimator may not be necessary. 

10 Alternatively, if the fitness measure indicates that R A does not sufficiently model the 
measured correlations, then the non-parametric estimator may be necessary to 
address the color that is not modeled by the parametric estimator. 

The above describes a method and apparatus for deriving an impairment 
correlation matrix R from first and second impairment correlation matrices to 

15 accurately track variations in modeled and non-modeled interference, such as self- 
interference, multi-user access interference, etc., in a wireless receiver. The wireless 
receiver may be disposed in any wireless communication terminal, such as a base 
station and/or a mobile station. A parameter processor 1 10 of the wireless receiver 
may use the derived impairment correlation matrix R to calculate various signal 

20 processing parameters. For example, using the final impairment correlation matrix R 
derived according to any of the above-described processes, weighting factors w* may 
be calculated in weighting calculator 118 according to any known method. For 
example, a weighting vector w comprising one or more weighting factors w k may be 
calculated according to: 

25 w = R l c. (Eq.7) 

Alternatively, the weighting vector w may be calculated according to: 

w = (f" 1 ^ Rd 1 F~ 1 c+ R a 1 c . (Eq. 8) 
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According to yet another method, the despread values are whitened and then 
combined using: 

w = Rd 1 F" 1 c. (Eq. 9) 

In any event, the weighting factors w* in the weighting vector w are used in the G- 
5 RAKE receiver to improve the performance of the receiver, as described above. 

Further, using the final impairment correlation matrix R derived according to 
any of the above-described processes, a signal-to-interference ratio (SIR) may be 
estimated by SIR calculator 1 19 according to any known method. For example, the 
SIR may be estimated according to: 

10 SIR est =c H R" 1 c (Eq. 10) 

SIR estimation is a part of the power control, which controls the transmit power 
necessary to maintain the communication link. SIR estimation is also used in link 
adaptation to provide the highest possible data rate for a given power budget. In any 
event, the final impairment correlation matrix R of the present invention may be used 

15 to improve the power control process associated with the transmitter by providing 
improved SIR estimates. 

The above describes an impairment correlation matrix estimation process and 
apparatus. To facilitate the above discussions, specific examples are provided. 
However, these examples are not intended to be limiting. For example, while a G- 

20 RAKE receiver is used herein to describe the invention, those skilled in the art will 
appreciate that the present invention is not limited to G-RAKE receivers, and is 
applicable to any spread spectrum receiver that estimates impairment correlation 
values. For example, the present invention is also applicable to the receiver 
described in U.S. Patent Application No. 10/672127 entitled "Method and Apparatus 

25 for RAKE Receiver Combining Weight Generation," filed 26 September 2003. 
Further, while the above described impairment correlation estimators 120 use 
parametric and non-parametric estimators, it will be appreciated by those skilled in 
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the art that any impairment correlation estimator may be used in place of either the 
parametric and/or the non-parametric estimator. The two impairment correlation 
estimators may be different types of impairment correlation estimators, as shown and 
described above, or they may be equivalent impairment correlation estimators 
5 programmed with different settings. Further, it will be appreciated by those skilled in 
the art that more than two impairment correlation estimators may be used according 
to the present invention to derive the final impairment correlation matrix R. 

In addition, those skilled in the art will appreciate that the above described 
impairment correlation matrix estimation process may be applied to wireless 

10 communication systems with one or more transmit and/or receive antennas. For 
example, a wireless communication system that includes multiple transmit antennas 
may use a receiver with an impairment correlation estimator 120 that includes a 
parametric estimator 1 30 for each of the transmit antennas. A wireless 
communication system that includes multiple receive antennas may have a non- 

1 5 parametric estimator system 1 50 that estimates impairment correlations across all 
fingers of the G-RAKE receiver, where each finger may differ in delay, antenna, or 
both. 

The present invention may, of course, be carried out in other ways than those 
specifically set forth herein without departing from essential characteristics of the 

20 invention. The present embodiments are to be considered in all respects as 
illustrative and not restrictive, and all changes coming within the meaning and 
equivalency range of the appended claims are intended to be embraced therein. 
Further, it should be emphasized that the term "comprises/comprising" when used in 
this specification is taken to specify the presence of stated features, steps, or 

25 components, but does not preclude the presence or addition of one or more other 
features, steps, components, or groups thereof. 
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